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SUMMARY 

An analysis was made of flight test data to obtain information 
representing the atmospheric turbulence Inputs to the test aircraft. The 
basis of the method was the use of the linearized equations of motion of 
the aircraft, with the addition of certain terms representing the aerody- 
namic effects of the turbulence on the aircraft. 

These equations were set up on an analog computer and verified 
by the technique of matching transient response flight data. Then, through 
the use of a feedback system, the desired turbulence quantities were ob- 
tained as system outputs. Inputs to the system were the aircraft respo/iso 
and control deflections, obtained from the flight test data. 

The method was found to be valid, assuming the aircraft analogue 
was correct. The resulting data appeared suitable for further analysis to 
determine its statistical qualities. 
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NOMENCLATURE 



PA Power approach 

Cl Lift coefficient, I i ft/qS 

Cm Pitching moment coefficient, pitching moment/qS 
C_0 Rolling moment coefficient, rolling moment/qSb 
Cn Yawing moment coefficient, yawing moment/qSb 

Cy Lateral force coefficient, lateral force/qS 
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1^ moment of inertia about x axis, slug - ft® 

ly moment of inertia about y axis, slug - ft® 

1 2 moment of inertia about z axis, slug - ft® 

S wing area, ft® 

V true airspeed, ft/ second 

Vq initial airspeed, ft/ second 

W aircraft weight 

H pitching moment of inertia parameter = 

mjxc® 

b wing span, ft. 

c wing chord, ft. 

c mean aerodynamic chord, ft. 

g acceleration due to gravity, ft/sec® 

h aircraft altitude, ft. 

m aircraft mass, slugs 

n normal acceleration, az/g 

3z normal acceleration, ft/sec® 

ay lateral acceleration, ft/sec® 

p = cp = rolling velocity, radlans/sec. 

q = Q = pitching velocity, radians/sec. 

r = i]f = yawing velocity, radians/ sec. 

t = time, seconds 

u = longitudinal velocity perturbation 

V = lateral velocity perturbation 

\'i = normal velocity perturbation 



I INTRODUCTION 



In accordance with a proposal submitted to the Bureau of Naval 
Weapons, the James Forrestal Research Center at Princeton University in 
late 1962 undertook an investigation into the lateral flying qualities of 
operational Navy aircraft In the PA configuration approach to carrier 
landi ngs. 

As noted in Reference I, the official proposal, it is typical 
that during this type of approach the lateral stability and control charac- 
teristics of some high performance aircraft deteriorate with increasing 
angle of attack and decreasing airspeed. This is common in many current 
sweptwing aircraft. The situation is aggravated by the added difficulty 
of approach and landing on aircraft carriers, where the pilot is required 
to maintain a flightpath with extremely narrow tolerances. The difficulty 
is further compounded by natural atmospheric disturbances, and those caused 
by the airflow about the carrier hull and superstructure, which are often 
encountered In the final few seconds prior to touch-down. 

It was therefore proposed that a flying qualities research pro- 
gram be conducted at Princeton, utilizing a variable stability aircraft to 
obtain pilot opinion ratings. The approach to a carrier was to be simula- 
ted as closely as possible, including effects of the forementioned atmos- 
pheric disturbances. These turbulence effects were to take the form of 
recorded signals applied as inputs to the electronic servos driving the 
control surfaces of the variable stability aircraft. It was therefore con- 
sidered desirable to represent the turbulence as accurately as possible 
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since this factor would probably have an important effect on pilot rating 
in this particular investigation. 

To this end, a suitably instrumented Navy fighter aircraft from 
the Naval Air Test Center, Patuxent River, Maryland, made a series of ac- 
tual carrier approaches under varying conditions of wind velocity across 
the deck and cross wind component. Recordings were made of the aircraft 
response in five degrees of freedom, together with recordings of the vai — 
ious control surface deflections. These data were analyzed on an analog 
computer setup of the aircraft equations of motion. By a process to be 
described, it was possible to determine the atmospheric disturbance infoi — 
mation In a form suitable for application to the variable stability aircraft. 
This report is concerned with the first phase of this research 
project: the acquisition, reduction, and analysis of the flight test data 

associated v/ith the carrier approaches made by the Naval Air Test Center 
ai rcraft. 
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1 1 EQUI PMENT AND PROCEDURE 
I) FLIGHT TEST EQUI PNENT 

The data acquisition flights were flown in an F4B fighter air- 
craft piloted by Navy test pilots from the Naval Air Test Center. The 
F4B (Phantom II) is a two-place, swept-wing, high-performance interceptor 
manufactured by the McDonnell Aircraft Corporation of St, Louis. It is 
powered by two General Electric J-79 turbojet engines, each producing a 
maximum of 17,000 pounds thrust at sea level. 

The aircraft was completely instrumented for the data acquisition 
by the Instrumentation Branch of the Flight Test Division, Naval Air Test 
Center, The rear cockpit seat was removed from the aircraft to provide 
space for the instrumentation package, A summary of the flight instrumen- 
tation follows: 

Power Sudd I V ! The basic aircraft power package was utilized as 
the instrumentation power source. This consisted of the normal 28 volt 
DC supply, and 115 volt, 400 cycle AC source. 

Sensing and Transducing : The physical quantities measured by the 

instrumentation system are shown in Table I, Potentiometer type transducers 
were chosen in order to provide compatibility with the telemetering system. 
The potentiometers provided signals of adequate magnitude, and through 
their excitation source also provided zero and full scale reference vol“ 
tage for the telemetering unit. 

Table I also lists the basic instrument used In the measurement 
of each quantity. The rate gyros which sensed the various angular velo- 
cities, and the potentiometers which measured the required position angles. 
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were all types commonly used in flight test instrumentation. It should 
be noted that the actual control surface deflections were measured rather 
than simple stick and pedal deflections, due to the non-linear character- 
istics of the stability augmentation system. Use of the actual deflec- 
tions permitted the aircraft equations of motion to be written in simple 
form. 

Airspeed and altitude were measured by W, 0. Leonard force balance 
type transducers. These units were located on the lower shelf of the rear 
cockpit instrumentation sled, as shown in Reference 2. The transducers, 
operating on the force balance principle, provided linear outputs for aii 
speed and altitude without the use of cams, function gears, shunts, etc. 

The force sensors were isolated from the output devices, thus minimizing 
response to transients. As a result, the units were characterized by a 
pressure accuracy of up to ,\%. In the case of the altimeter this meant 
accuracy to within three inches for a two hundred foot aircraft altitude. 

Lateral and normal acceleration were sensed by Donner model 4310 
linear accelerometers installed with a range of + .5 G and + 2.0 G, respec- 
tively. These units exhibited natural frequencies much higher than the 
test aircraft, thus precluding resonance problems. Their response was es- 
sentially flat to 8 cps. It was determined, however^ during installation, 
that some H. F. system noise was being passed by the accelerometer. In 
order to eliminate this noise, a low pass filter was therefore designed 
and installed between the accelerometer output and the telemetering unit. 

The filter, designed for a 5 cps corner frequency, eliminated 
the high frequency noise but provided essentially flat response over the 
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expected aircraft response range, which was estimated to be no higher than 
I cps, well above the Dutch roll frequency. 

A schematic of the filter circuit is shown below, together with 
a sketch of the filter frequency response characteristic: 
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Telemetering and Record! no ; Telemetering equipment of the pulse 
width type was installed in the aircraft as part of the instrumentation 
package. This equipment, manufactured by the Applied Science Corporation 
of Princeton (ASCOP), was the manufacturer's D Series Pulse Width Multi- 
codei ^Telemeter, capable of sampling up to 43 channels at a sampling rate 
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of 20 times per channel per second, using a Model 0SG20 Commutator-Gear 
Box. 

Since each channel was sampled 20 times per second, it was felt 
that good fidelity would be obtained with input frequencies as high as 
the 2 to 5 cps range. Many of the measured quantities were sampled on each 
of two or three channels, and in these cases it was felt that good repro- 
duction was assured at all frequencies encountered in the investigation. 
Table II lists the measured quantities and the number of sampling channels 
used for each. For those quantities sampled in more than one channel, the 
sampling was, as far as possible, equally spaced in time. 

The sampled data were then converted into RrJ data by the Model 
DKH-I Keyer. 

The pulse width information was then recorded directly in the 
aircraft by means of an Ampex Model 800 tape transport which was mounted, 
together with its associated power supply, in the rear cockpit installation. 

Provisions were also made for simultaneous voice recording by 
the test pi lot. 

2) GROUND EQUIPMENT 

The following units were available at the Forrestal Research 
Center for analysis of the flight test data: 

Magnetic Tape Recorder : An Ampex Model 309C dual track, three 

speed tape recorder was available for playback of the recorded flight data, 
and also for use in the preliminary analysis studies which were conducted. 

Telemetering Ground Station : The telemetering unit used for de- 

coding the taped pulse width flight test data was the ASCOP M Series PW 
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Ground Station. This unit continuously sampled the recorded pulse width 
information and converted this data into continuous voltage form. The 
information in this form was then compatible with the analog computer for 
analysis purposes. 

GEDA Analog Computer ; The basic tool in the analysis process 
was the Model L3 (GEDA) linear electronic differential analyzer, manufac- 
tured by the Goodyear Aircraft Corporation. This computer provided twenty- 
four automatically stabilized DC computing amplifiers, each with an open 
loop gain of greater than 5 x lO’’, and with negligible drift. The computer 
incorporated an automatic error indicator, and had a guaranteed accuracy 
of 1^. There were provisions for accurately setting computer board po- 
tentiometers using a special null indicator. 

Data Recorder ? Visual presentation of data was made possible 
by a Sanborn Model 154-lOOB four channel recorder. This instrument had 
a very high natural frequency (42 cps), and low frequency response flat 
to zero cps. 

3) INSTRUMENT CALIBRATION 

Calibration of the flight test instrumentation was conducted by 
the Instrumentation Branch of the Flight Test Division, Naval Air Test 
Center. The calibration charts are shown in Reference 2. All gradients 
were determined to be essentially linear within the range of measured data 
encountered during the test flights. 
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4) PROCEDURE 
Data Acquisition ! 

In December, 1962, the instrumented F4B aircraft, piloted by a 
Navy test pilot, made a series of thirty-three landing approaches to the 
aircraft carrier USS Enterprise. These approaches and landings were made 
over a span of five separate flights. The heading and speed of the car- 
rier were varied during the course of the series to provide a number of 
different conditions of wind velocity across the deck, and cross wind com- 
ponent. The wind conditions for each approach are shown in Table III. 

It was expected that the most severe turbulence conditions would 
be experienced in the presence of a "right to left" crosswind. This was 
due to the extremely turbulent flow about the carrier "island" being direc- 
ted into the path of an approaching aircraft, usually at the critical point 
occurring a few seconds prior to touchdown. It was felt that by making 
the approaches under varying wind conditions, it would be possible to ob- 
tain turbulence data which would be representati ve of the full range of 
severity normally encountered. As noted in Reference 3, the turbulence 
"burble" associated with the Enterprise was probably the extreme, due to 
the very box I ike design of the ship's "island". 

The Naval Air Test Center provided facilities aboard the Enter 
prise for tape playback after each test flight. On the spot assessment of 
data quality was thus possible, and minor adjustments were made, as neces- 
sary, in the instrumentation. 
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Aircraft Calibration : 

Since the basis of a successful analysis of the flight test data 
was to be the setup of the aircraft equations of motion on the analog com- 
puter, it was necessary, as the next step, to obtain a valid set of equa- 
tions representing the F48 aircraft in the PA conf iguration. The basic 
method used was the matching, on the computer, of transient response flight 
test data. The method will be discussed in some detail in a subsequent 
section. 

To this end, two additional flights were made of the test aii — 
craft for the purpose of obtaining transient response data. The aircraft 
was flown in the PA configuration, and pilot control inputs of the pulse, 
step, and doublet types were made. The aircraft response was recorded 
as before. 

In order to prevent stability augmentation equipment inputs to 
the system, these units were made inoperative during the transient res- 
ponse data runs. 

Data Ana I vsi s ; 

Having thus acquired the necessary flight test data, the actual 
analysis was then conducted at the Forrestal Research Center. This analysis 
is discussed in detail In the following section of this report. 
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in DATA ANALYSIS; OESCRfFTI ON AND DISCUSSION 

The first step in the data analysis was the determination of a 
valid set of differential equations representing the F48 aircraft in the 
PA configuration. It was assumed that the motions of the aircraft whi le 
approaching a carrier in this configuration could be described by the usual 
linear equations. This assumption seemed valid since the final part of 
a typical carrier approach, where the turbulence effects are normally most 
severe, is basically a constant power, constant airspeed, shallow descent 
to a touchdown; the aircraft motions are approximately perturbations from 
the steady state glide path. 

The resulting equations are shown below, written in dimensional 
form for compatibility with the flight data and analog computer. For con- 
venience, the terms of each equation have been divided by the coefficient 
of the highest appropriate derivative. The assumption was also made that 
the aircraft forward velocity remained essentially constant during the 
runs, thus eliminating the drag equation, and also eliminating longitudinal 
velocity terms from the remaining equations. 
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These uncoupled sets of equations were then separately mechan- 
ized on the GEDA analog computer. Schematic circuit diagrams of the two 
setups are shown in Figures 1 and Z. 

All factors in these equations were known accurately, except 
for the various stability derivatives. Approximate values for the deriva- 
tives, for the PA configuration, were obtained from a report of the 
McDonnell Aircraft Corporation. These values determined the initial po- 
tentiometer settings for the computer setups. Shown below is a I i st of 
the various quantities utilized in the initial formulation of the aircraft 
equations for the PA conf iguration: 
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134 knots = 224 feet per second 

34000 lbs. 

16 degrees 

1. 1 

.002378 slug/ ft® 

39 feet 
530 feet® 

16 feet 

23880 slug - ft® 

121007 slug - ft® 

131 177 slug - ft® 
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The equations of motion based on these values are as follows: 
Longitudinal i 

1. a - 9 + .406 a + .0436 63 = 0 

2. 0 + . 162 0 + .067 a + .357 a + .792 63 = 0 

3. n + 2.81 a - .31 1 6s = 0 

Lateral ; 

4. - p - .0832 p - .990 } + .0 109 9 + . 139 9 + .0166 6p = 0 

5. - 9 - I .25 9 + .986 i - 12.2 p - .510 8^ + 2.80 6g = 0 

6. - J - .264 ^ + .0288 9 + 1.89 p - .693 6^. - .0556 63 = 0 

Improving the validity of these equations as representations of 
the aircraft thus required a more accurate determination of the stabi llty 
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derivatives. This was accomplished by a transient response curve fitting 
technique, such as outlined in Reference 4, and described in detail in 
Reference 5, The mathematical basis for the validity of the method is 
found in Reference 6, 

For this method a number of the transient response runs were 
selected by reference to Sanborn recordings of the aircraft response. The 
selection was made on the basis of freedom from telemetering dropouts, and 
relative freedom from apparent extraneous inputs# such as turbulence effects. 
These response curves were of the rudder pulse, rudder doublet, ai leron 
step, aileron doublet, stabi lator pulse, and stabilator doublet types. 

Transparent overlays were made from the Sanborn records, incor- 
porating the scale changes necessary to provide agreement with the scaling 
of the computer outputs. It was then possible, by means of the ground 
station equipment, to apply voltages to the aircraft analogue represent- 
ing the actual aircraft control motions, as forcing functions to either the 
longitudinal or lateral sets, as appropriate. Sanborn recordings of the 
analog responses were made, all to the same scale as the transparent over- 
lays. The two sets of data were compared, and by an iteration procedure, 
the potentiometer settings representing the various stability derivatives 
were changed until close agreement was attained between flight response 
and analog response. 

In the lateral case the matched quantities were roll rate and 
yaw rate. It was found that a close match to the flight data could be 
fairly rapidly obtained by corrections to those derivatives which had major 
effects on the important lateral response parameters, i.e., Dutch roll 
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period and damping, initial roll and yaw magnitudes, spiral stability. 
Improvements to the initial match were obtained by a systematic iteration 
procedure, varying the potentiometer settings in turn until the desired 
match resulted. 

In the initial study of the lateral responses it was immediately 
apparent that the flight data was not ideal. During even the best runs, 
the aircraft obviously had not been In steady flight at the instant of the 
control inputs; there were initial conditions which had to be considered 
in the matching process. By means of the analog initial condition cii — 
cuitry it was possible to apply to the computer, as Initial conditions, 
the actual flight data response in roll rate and yaw rate existing at the 
instant of control deflection. This was considered to be an improvement 
but it was noted that since the lateral equations represented a fourth 
order system, two additional initial conditions, P and cp, for example, 
were necessary for accurate solution. These were not aval lable and some 
error was thus to be expected. In the iteration procedure it was found 
that the best matches were obtained with flight data where the aircraft 
was very close to steady state flight conditions at the time of control 
deflection. This was to be reasonably expected since in these cases all 
initial conditions were approximately zero. 

It was also apparent from observation of the flight data that 
some atmospheric turbulence existed during the tests. This, of course, 
could not be accounted for in the matching process, and was thus a source 
of error. For this reason the smoothest runs were selected for the match- 
ing process. 
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It was noted also that during the various transient response runs 
there were inadvertent inputs of those control surfaces which theoretically 
were being held in the neutral position* These inputs were also played into 
the analog. 



It can be seen from Table I that only the left hand ai leron and 
spoiler were instrumented for deflection measurement. Since the equations 
of motion as formulated on the computer were written in terms of aileron 
inputs, which are the sum of right and left ai leron deflection, it was 
necessary to derive some function of left spoiler deflection which might 
be utilized as equivalent right aileron deflection. To this end, a cali- 
bration run was made to determine the relationship between ai leron and spoiler 
deflection. The results were of the form shown in Figure 3. The relation- 
ship between the two was shown to be approximately linear. Therefore, it 
was possible to express right aileron as a linear function of the recorded 
spoiler deflection. The aileron signal fed to the analog was; 



(7) 



6;, = 6, + = 6^ - K 6 






It can be seen from Figure 3 that each aileron exhibited a small 
upward deflection range of approximately three degrees. This property was 
ignored in determining the proportional i ty constant between ai leron and 
spoiler. It was felt that the upward deflections were primari ly for the 
purpose of assisting the effect of the spoi ler deflections on the same 
wing. In the slightly open position the spoilers were assumed to be rela- 
tively ineffective in disrupting the airflow, and the small upward aileron 
deflection made up for the deficiency, thus maintaining system linearity. 
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By means of the analog limiter circuits, the voltages representing these 
upward aileron deflections were eliminated from the computer inputs. 

It should be noted that in spite of the turbulence and initial 
condition problems. It was possible to obtain good matches between flight 
data and analog output. The accurate determination of all the various 
stability derivatives would probably have been difficult because of these 
problems, but the primary interest here was the matching of aircraft re- 
sponse, This could be done by matching the important lateral parameters 
mentioned previously. Since these parameters are primari ly controlled by 
relatively few of the derivatives, it was possible to obtain a good match 
even though some of the derivatives might have been in error. 

It should also be noted that since the instrumentation used In 
obtaining data for this matching was the same as subsequently used during 
the actual carrier approaches, then to the extent that any i nstrumentation 
errors would have caused false matching, the effect would tend to be can- 
celled during the actual approach data acquisition. 

Shown below are the circuits utilized to provide inputs to the 
analog from the recorded data. By means of appropriate potentiometer set- 
tings the quantities were all scaled correctly for compatibility with the 
analog setup. Biasing of all quantities to the correct zero position was 
done by adjusting the biasing knobs of the ground station while referring 
to voltmeter readings of the various tape outputs. 

Shown in Figures 4, 5, and 6 are some of the matched transient 
response curves. It i s to be noted that the aileron step inputs were used 
primarily to match initial roll rate amplitude. It was felt that the aircraft 
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was out of the linear range shortly after application of an aileron step 
and this precluded any additional useful matching. 
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The longitudinal responses were matched by comparing pitch rate 
and normal acceleration. The method utilized and the problems encountered 
were similar to those discussed in reference to the lateral equations. 

The phugoid mode was not visibly excited in any of the transient runs and, 
therefore, the matching was confined solely to the short period oscil- 
lations. The matching iteration procedure converged much more rapidly than 
in the lateral case. This was to be expected since in this longitudinal 
approximation there were only two degrees of freedom and, consequently, 
fewer controlling derivatives to be adjusted. 
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Shown below are the input circuits to the longitudinal analog; 



ffu/h. t 

- 



3Cd Uno 

"Qo:^ 



~ ^CO- 

J 



o 



/ 

-AAAA 







y- 



4 



fo IC // 



In Figures 1, 8, and 9 are shown some of the matched transient 
response curves. 

It should be noted that the initial matching, both longitudinally 
and laterally, showed a number of di screpancies between the actual aircraft 
response and the analog response which was based on the stability deriva- 
tives obtained from the McDonnell report. In the lateral case the analog 
response showed a predicted Dutch roll period of approximately 4.5 seconds, 
with neutral damping. As can be seen in Figures 4 and 5, the aircraft 
actually exhibited a well damped Dutch roll of about a 3 second period. 
Longitudinally, the major discrepancy was In the apparent degree of aii — 
craft static stability, which was considerably greater than the predicted 
value. 

These discrepancies are illustrated more clearly by reference to 
the final equations of motion which resulted from the matching process. If 
these equations, which are shown below, are compared term by term with the 
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original set, the changes which had to be made during the matching are 
apparent. 

Longitudinal t 

( 8 ) a - 0 + .472 a + .0436 65 = 0 

(9) 0 + .838 0 + .067 a + 9.00 a t 4.50 6 g = 0 

(10) n + 3.25 a - .311 6 s = 0 

Latera I ; 

(11) - 8 - .0956 8 - .990 Jr + .0110 cp + .140 cp + .0209 6 p = 0 

(12) - (p - 2.48 - .984 \ - 15.8 8 + .497 6 p t 5.48 6 g = 0 

(13) - \jr - .569 i + .0288 9 + 3.38 8 “ .674 6 r - .061 6 a = 0 

It can be seen that in the lateral set the directional stability 
and yaw damping, which are the primary parameters affecting Dutch roll 
period and damping, respectively, were each approximately doubled in the 
matching process. Longitudinally, it can be seen that the static stability 
was also increased considerably during the matching. Although the high 
value of Ctn^, which finally produced a good match, was probably unrealistic 
it at least indicated that the predicted value was too low. It is noted 
again that accurate determination of each stability derivative was by no 
means assured. However, the aircraft response curves were fairly accurately 
matched, and this was the important point for purposes of this Investigation. 

With the equations of motion of the aircraft thus formulated, 
the next step was the actual analysis of the flight data to obtain the 
turbulence information. Preliminary studies were made to determine the 
feasibility of the proposed method of solution. For these studies the 
lateral equations of the Navi on aircraft, as obtained from Reference 5, 
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were set up on the analog computer. The basis of the method was the well- 
known feedback principle: 




From the above diagram, C/R = A/(l + AB), If A, the forward 

loop gain, is very much greater than unity, then, 

C « 1 

R B 

If B were an aircraft transfer function, e.g., 9/63/ then, 

C, fti I 

R <^)/ 6 a 4 * 

It seemed possible, then, that if known aircraft response in roll rate, for 
example, were fed into such a circuit at "R", where "B" was the aircraft 
transfer function as set up on the analog computer, then the circuit out- 
put at "C" would be a voltage representing the aileron deflection which 
had caused the roll response. 

The validity of this theory was tested by recording on tape the 
analog response i n \jr, and p, to an arbitrary sequence of voltages repre- 
senting aileron inputs to the system. This was accomplished by connect- 
ing a potentiometer to the appropriate points of the analog circuit and, 
in effect, "flying" the aircraft in the lateral modes. The previously 
discussed feedback circuit was then connected as shown below: 
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The taped response in roll rate was then played into the circuit at "R" 
and the output "C" was compared with the original aileron signals. The 
match was excellent with forward loop gains as low as 10. Below this value 
the reproduction decreased in accuracy since (C/R) (I/B) was no longer 

a good approximation. 

The next logical step was to study the effects of two control 
inputs, both aileron and rudder. Again the analog was "flown", this time 
with aileron and rudder. The same aircraft response quantities were again 
taped. The analysis circuit was similar to the previous case, with the 
addition of another loop. The circuit Is shown below in block diagram 
form. The taped responses In 9 and ^ were played Into the circuit at R 
and R' respectively, and the outputs at C and C' were compared with the 
original 6 g and 6 p. Again the match was. excellent, down to forward loop 
gains of 10 . 
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Figure 10 shows a more detailed schematic of this last analysis 
circuit. It was important to provide negative feedback in the loops in 
order to maintain system stability. 

The results of this preliminary study pointed up the uniqueness 
of the solutions to the aircraft equations of motion. A given combination 
of aileron and rudder inputs results in a unique combination of roll/ yaW/ 
and sideslip response. 

Originally/ it had been proposed that this exact analysis method 
be employed in the study of the flight test data. In the lateral case/ 
for example/ aircraft response in roll rate and yaw rate were to be applied 
to a similar analysis circuit/ where the analog would be of the F4B, From 
the resulting outputs of 6p and 6g then would have been subtracted the 
recorded flight test control deflections. The difference would be an 
aileron and rudder signal lepresenting the turbulence input to the F48. 
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This information, suitably modified to account for the difference in con- 
trol power between F4B and Navion, might then be applied to the variable 
stability aircraft as a simulation of the turbulence input. 

The longitudinal analysis would be similar. Aircraft response 
in pitch rate would be applied to the appropriate circuit and the output 
6 g would yield an equivalent 65 representing the turbulence input. 

This method, however, obviously neglects the linear perturbations 
of aircraft motion due to turbulence, those motions along the aircraft 
axes. It assumes that the atmospheric turbulence effects may be des- 
cribed accurately by the resulting moments about the aircraft axes. For 
purposes of this research project reasonably valid simulation may be at- 
tained based on this concept. 

However, in the interest of more accurate gust simulation, it was 
decided to conduct some initial research into the possibility of obtaining 
turbulence data which might be expressed in five degrees of freedom, thus 
including any lateral and normal acceleration effects. In addition, the 
possibility existed that some qualitative information concerning the chai 
acter of the turbulence itself might result from such a study. 

It was known that it is possible to visualize the atmosphere 
through which the aircraft moves as being in motion at some steady mean 
value, with random turbulent fluctuations superimposed. These random velo- 
cities may be considered as inputs to the aircraft linear system, which 
may respond in six degrees of freedom. The velocity inputs are assumed 
small, permitting linearization of the resulting aerodynamic forces. 
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As outlined in Reference 1, it is convenient to introduce these 
inputs into the equations of motion in the form of velocity components 
Ug, Vg, and Wg (subscript "g" standing for gust). If it is assumed that 
the gust wavelengths are all large relative to the aircraft dimensions# 
then the gradient, or velocity variation, of these gusts across the length 
and span of the aircraft may be considered essentially linear. It can 
be shown that under this assumption the effects of the gust components and 
their gradients are aerodynamical ly equivalent to perturbations In angle 
of attack, sideslip, roll rate, yaw rate, and pitch rate. These are noted 
as ag, 3 g, cpg, fg, 0g. 

The aircraft equations of motion incorporating the modified var- 
iables thus become; 

Lonoi tudi na I ; 

• * c c 

( 14 ) a - 0 + Jict. (a + a„) + s 6s =0 
2 T ^ 2t ■ 



(15) 0 - (d + eg) - a - 22a (a + ttg) - ^"6 s 65 =0 

Ht ^ Ht ^ 

Latera I ; 

( 16 ) 3 - (3 + 3q) + ^ - (rjr + | ) + 2 a 9 “ tcp + cpg) 

2t ^ 4p, 2 t4|x 








2t 



2t 



0 
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(17) 9 ■ <9 + 9g) “ (t + iifg) ” (p + Pg) - 6 



2JxT 



2JxT 






JxT* 



_ P^i6p 
JxT® 



6 p = 0 



(18) J - £2I_ (t + L) - ^ (9 +9^) - (P + Pc) - !fn^ «e 



2 J 2 T 



2 JzT 






6p = 0 



With the gust terms rearranged to point out their roJe as forcing functions, 
the equations are written as: 

Lonoi tudi na I : 

(19) d “ 6 + a + **^5 63+ ^*^0. aq =0 
2 T 2 t 2 t 



( 20 ) 9 - d - a - 65 ~ §c " 2 l]a aq =0 



Ht 



Ht 



H-r 



Ht® 



Ht 



Ht* 



Lateral : 



( 21 ) P - ^ P + i - 2 IL ^ + 2 i 9 - ^ 9 - 6 a - 6 r - 2 ^ P 



2 T 



4 u. 



2t 



- "'Vp 9 - 


2^63 


4 Ui 


2t 


^Yr ; 


^Yp • 
^9g 


4 p 


4 iJ, 


_ PC^6a 6 


- 



2t 



= 0 



2 T 



2 JxT 



2 JxT 



JvT* 
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XF’ 



2 JxT 



=0 
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(23) 



^ cp - p - 6a - 



2J2T 



2JzT 



JzT^ 



JzT" 



J7T® 



6r - ^nr 
2J2T 



2J2T ^ J^T® 



= 0 



It i s to be noted that these equations are not valid for short 
wavelength gust components. 

However, it can be shown that the effect on aircraft response of 
the higher frequency gust components is negligible. For example, the ef- 
fective gust frequency, Og, may be expressed by; 




where \ is the gust wavelength. For a short X, say equivalent to the aii — 
craft wingspan, Wg is on the order of 5 cps, for the F4B in the PA configu- 
ration, This is very much higher than the aircraft natural osci I latory 
frequencies and consequently the aircraft response is negligible. 

In addition, for gust wavelengths of this order, and shorter, 
the aircraft begins to "average out" the gust velocities, with the result- 
ing net gust acting on the airframe being near zero. Thus, for these higher 
frequency components there is a compounded attenuation - that due to the 
aircraft dynamic characteristics, and that due to the cancelling effect 
just described. 








liii . 
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|t should be noted that If the F4B were to respond dynamically 
to the higher frequencies, this last effect, that of cancelling, would 
cause some error if the acquired turbulence information were applied to the 
study of an aircraft differing greatly in size from the F4B, It is quite 
possible that the aircraft in question might respond to those gusts which 
were cancelled by the F48 and were thus missing from the acquired data. 
However, all aircraft which will be simulated by the variable stability 
Navion, for this investigation, are similar in size to the F46, and there- 
fore could be expected to "average" the gusts in the same manner. 

With these modified equations of motion in mind, it was desired 
to obtain, as a solution, the various gust quantities. The method, although 
based on the previously discussed feedback principle, can probably be best 
described by a single degree of freedom example. For this purpose it is 
assumed that the aircraft can respond only In roll, and that the only gust 
quantity present is <J)g. The equation of this system is then; 

(24) 9 - ^^P 4 , - 6 a “ cp = 0 

2J^t JxT® 2J^r ^ 

This equation is set up on the analog computer and a feedback 
circuit is connected as shown on page 29. 

The output of the forward loop amplifier is defined as - cpg, 
and this quantity, modified by the appropriate proportiona I ity constant, 
is used as an analog forcing function. The recorded values of - and 
- 6a are also used to force the computer setup. The following equation, 
which expresses the aircraft response, is satisfied at point A: 

"^recorded ~ 'P control response ■*' “Pgust response 
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Assuming that the circuit is dynamically stable, then the - cpg 
quantity actually appears at the amplifier output. This circuitry is simi- 
lar to those previously discussed. 

The more complicated circuit which was proposed to solve the en- 
tire lateral set is shown in detail in Figure II. It is, however, based 
on the same ideas. The three equations which are satisfied at the forward 
loop amplifier inputs are; 

^recorded *" ** ^gust response *" ^control response ” ^ 

• • • —A 

'Precorded ” 'Pgust response “ 'Pcontrol response 

• . , 

Vecorded “ 'I'gust response ^control response “ ® 

It is to be noted that in the case of 3, the recorded quantity 

from the tape was based on an accelerometer measurement and, thus, con- 
tained the quantity 3 g, in addition to response^ aircraft re- 

sponse to 3 g, Therefore, to satisfy the above equation, the Pg was fed 
back to the amplifier input. In the case of 9 and the aircraft response 
only was included in tho measui ed response since the sensors were rate 
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gyros, Approximately i./>finite resistance feedback was therefore provided 
around the forward loop amplifiers in these two channels. 

It was discovered, however, that the circuit shown in Fi gure 
II was not satisfactory. During operation of the circuit' the output quan- 
tity - Pg exhibited a di vergence„. This Indica-ted some positive feedback 
within the multi loop systenv It can be seen that There were nine loops 
through the aircraft analog* IT seemed probable that the - Pg outpvt was 
controlled primarily by the loop through amplifier 9 which is shown below 
schematical ly: 




This circuit obviously provides positive feedback and is thus 
unstable. This problem emphasized the point that although the equations 
were physically satisfied in the circuit as shown, it was necessary to 
insure system stability in the electronic sense. A study of the circuit 
showed that there were additional unstable loops, but that their effect 
was minor compared to the effect of the stable loops which fed back to 
the same points. 

In order to provide stability, the analysis circuit was modified 
to the form shown in Figure 12. The difference here is that the - Pg quan- 
tity input to the analog is now taken directly from the recorded data. The 
analog circuit was modified to maintain satisfaction of the equations of 
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motion. Circuit stability was now attained since the unstable - 3g loop 
through amplifier 9 was now eliminated. 

To test the validity of the analyzer circuit, a preliminary study 
was made in which the aircraft analogue was forced by a ficticious 3g quan- 
tity as obtained from a potentiometer output. The aircraft response in 
roll rate, yaw rate, and sideslip was recorded on tape. It is to be noted 
that the recorded sideslip was made the sum of the analog 3 response and 
the forcing function, 3g. This was to provide similarity to the actual 
flight data where, as previously mentioned, the 3 response was based on 
an accelerometer measurement. These recorded responses were then played 
into the analyzer circuit in an attempt to duplicate the original forcing 
function. 

Figure 13 shows the original forcing function and the recorded 
aircraft responses. When these responses were played back into the ana- 
lyzer circuit, the cpg and ijig outputs were essentially zero, as expected. 

The 3g output was found to be a good reproduction of the Input wave form, 
although there was a very low frequency oscillation superimposed on the 
trace. This is shown in Figure 14(a), where the solid line Is a tracing 
of the original 3g input, |t was determined experimentally that this low 
frequency drift was strongly affected by the biasing adjustment of the 
tape output from the ground station equipment. Careful adjustment of the 
bias practically eliminated the drift; but since this method would not be 
feasible with actual flight test data, another method was sought. It was 
decided to utilize a high pass filter to eliminate the low frequency os- 
cillation. This was connected as shown: 
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Since the use of such a filter would result also in the attenua- 
tion of some of the higher frequencies, the corner frequency was selected 
on the basis of compromise between minimum attenuation and minimum divei — 
gence. The results are shown in Figure 14 (b), (c), and (d), A corner 
frequency of 1/44 cps, as provided by Rf = 7 meg Q, seemed to give the 
best reproduction. 

Throughout this discussion it has been assumed that the aircraft 

« 

response data were available In the form of 3, 9 , and It can be seen 
from Table I that lateral acceleration, ay, was recorded during the flight 
tests, rather than sideslip. It was necessary to obtain 0 from this data 
by means of the following form of the side force equation; 

(25) + p = m I ay-^^q>-2XrJ^ir - 6^ 

(p/2)Vo^S Cyp Cyp 2Vo Cyp 2 Vq Cyp 

The summing and scaling circuit utilized to solve this equation 
for application to the analyzer is shown below; 
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Si nee an accurate representation of p depended on reliable valuer ■ 
of Cyp, since these derivatives were not involved 

in the previous matching process, it was necessary to check their validity. 
The method used was to compare the quantity ay from the flight data with 
the derived ay obtained by summing the appropriate analog outputs in re- 
sponse to a rudder pulse. The results are shown in Figures 4, 5, and 6. 

It was felt that the matching was accurate enough to justify the use of 
P, as obtained from equation 25. The resulting equation for ay is: 



(26) + ay = (p/2)Vn^S Cw B+ (o/2)Vn^S b C„ (P/2 )Vq^S b Cy i 

P m 2Vn P 



m 



m 2Vn 



+ (P/2)^o S Cy 6r 
m ^ 
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The summing and scaling circuit for obtaining ay for the matching 
process is shown below; 
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The analyzer circuit designed for use with the longitudinal equa- 
tions was based on the same principles which governed the lateral circuit. 
The schematic of this circuit is shown in F'gure 15, The following equa- 
tions were satisfied at the inputs to the forward loop amplifiers: 

• . • 

^aircraft response “ ®gust response ^control response 

^aircraft response ~ *^gust ®gust response ^control response 

It should be noted that the reliability of the analyzer output 
was prpbably improved if the aircraft maintained a very steady, smooth 
flight path during the carrier approaches. The reason for this is that 
with the aircraft responses very small, the aircraft analogue, in the ana- 
lyzer feedback toops^ was relatively inactive. This was beneficial since 
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the analogue had known inaccuracies. It was therefore desirable that the 
aircraft stability augmentation equipment be operative as an aid in main- 
taining steady flight during the carrier approaches. 
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IV RESULTS 

The next step in the data analysis was to play Into the analyzer 
circuits the appropriate quantities from the flight test tape recordings. 

Of the five flight test tapes, one was determined to be virtually unusable 
due to superimposed noise. The four remaining tapes were each played Into 
the analyzer and Sanborn recordings were made of the five gust quantities 
for each carrier approach. Also recorded for reference purposes were the 
aircraft response in five degrees of freedom and the various control 
def lections. 

Figures 16 through 31 are reproductions of gust quantities from 
the lateral and longitudinal analyzers. Also shown are traces of aircraft 
response and control deflection. Aircraft altitude also appears on the re- 
cordings as an aid in orientation. It was apparent from the traces that 
some drift was still present in the pg channel, probably due to incorrect 
biasing at the ground station output. Further increase of the filter coi 
ner frequency was not practical due to the extensive attenuation of the 
higher frequency components. 

Those approaches which are missing from the set of Sanborn re- 
cordings were unusable due to the effects of "dropouts" during the analyzer 
operation. As mentioned previously, the tape of the second flight, con- 
taining approaches 8 through 13, was completely unusable due to poor quality 
of recorded information. 

Study of the traces revealed a number of points. It was expec- 



ted that if the analyzer circuits were completely valid the output gust 
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quantities would be completely random in character, reflecting the nature 
of the actual turbulence. However, Figure 19, for example, shows that 
there were occasional sections of the various gust traces which were perio- 
dic in character. Conversely, many of the remaining traces. Figure 16 for 
example, appear quite random. Some of the oscillations appearing were of 
the approximate character of the aircraft short period and Dutch roil oc- 
ci Nation, The implication here was that the F4B analogue was not com- 
pletely accurate. As a result, some of the aircraft response characteris- 
tics appeared in the analyzer output. 

In spite of the apparent imperfections in the turbulence data, 
some qualitative information was immediately available. The records showed 
that for all approaches there was a marked increase in the magnitude of the 
gust quantities in the few seconds prior to aircraft touchdown. This was 
a verification of the existence of the turbulent "burble" immediately astern 
of an aircraft carrier deck. 

In order to detect any additional trends In the turbulence quan- 
itites, the RMS value of each quantity for each approach was determined 
with the aid of the following circuit, which included the Donner Model 
3751 Function Generator: 
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(27) RMS 




/X(cin^)dt 



In order to eliminate the effect of any DC bias or low frequency 
oscillations, the various gust quantities were channeled through a high 
pass filter prior to entry into the above circuit, A corner frequency of 
.20 cps was selected for the filter, since it was felt that frequencies 
below this would probably not have a significant effect on pilot opinion. 

The RMS values of the various quantities were tabulated and are 
shown in Table II. Omissions in the Table were due to the effects of tele- 
metering dropouts. These tabulated values do not include the effect of the 
final "burble". No apparent trends were determined from study of the RMS 
tabulation, except for an indication of the general level of each quantity. 
These levels remained approximately constant during the course of the ap- 
proaches, indicating that atmospheric conditions also generally remained 
the same. 

It was also noted, with reference to Table III that the angular 
velocity gust quantities, ^g, ij»g, and 6g were all of approximately the 
same magnitude. The angular quantities pg and ag were also approximate ly 
equal. This was a strong indication that the turbulence away from the 
Immediate vicinity of the ship was indeed isotropic, as is normally assumed. 

It was expected that the magnitude of the gust quantities within 
the "burble" would be a reflection of the existing relative wind conditions 
over the carrier deck. Although no formal tabulation was made of the RMS 
values of the gust quantities within the "burble", study of the traces 
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revealed that these RMS values were approximately doubled in the final 
three seconds of each approach. There was no apparent difference in the 
character of the "burbles" with varying relative wind conditions. 

It was noted, however, from study of the gust traces that the 
"burble" was initially encountered at varying distances from the ship. 

Figure 32 shows the point of entering the "burble" for all of the approaches 
represented by Figures 16 through 31. It is interesting to note that those 
points furthest from the carrier were in the presence of high relative 
wind, with no x-wind component. Those points nearer the ship were with 
less wind velocity. Finally, those points appearing 2 and 3 seconds prior 
to touchdown were in the presence of "right to left" x-wind, indicating 
that the turbulence trail was being blown across the flightpath at this 
point. 

In Figure 33 are plotted the pertinent lateral mode data from a 
typical approach. The increase in the magnitude of the gust quantities in 
the final few seconds before touchdown is apparent. However, it was also 
noted that the gust quantities at this point had some periodic characteris- 
tics. In addition, the oscillations are approximately of the same period 
as the aircraft roll and yaw response. In fact, in the case of and ^g, 
the oscillations are exactly 180° out of phase. Because of apparent "im- 
purities" such as this in the turbulence data It was felt that a more de- 
tailed statistical analysis of the data should be made to determine its 
true spectral character. Such an analysis is beyond the present scope 
of this project. 
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V CONCLUSIONS AND RECOMMENDAT I ONS 

Based on the results of this investigation, it is concluded that; 

1. Accurate representations of the atmospheric turbulence may be obtained 
by the analysis method described in the report, assuming that the aircraft 
analogue is valid. 

2. The consistent simi larity in the magnitude of the angular velocity 
gust quantities, and the simi larity in the magnitudes of gust angular quan- 
tities, as shown in Table IV, indicated that the turbulence away from the 
immediate vicinity of the carrier was isotropic in character. 

3. The atmospheric turbulence in the vicinity of an aircraft carrier is 
character ized by a strong "burble" in the vicinity of the approach end of 
the carrier deck. It is not possible, based on these results, to give 

a quantitative description of this turbulence although it appeared that t/ie 
character of the "burble"did not vary noticeably with changing relative 
wind conditions. However, the point at which the "burble" is initially 
encountered is a function of the relative wind conditions. 

It is recommended that: 

1. A detailed spectrum analysis be made of the turbulence quantities re- 
sulting from this investigation, in order that undesirable portions of the 
spectrum be eliminated, making the information suitable for application 

to the variable stability Navion. 

2. A further statistical analysis of this data be made to determine quan- 
titatively the character of its components. 

3. An investigation be conducted into the applicability of this analysis cii 
cuit for the general solution of complex inverse transfer function problems. 
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TABLE 


1 


RECORDED 


FLIGHT TEST 


DATA QUANTITIES 


QUANTITY 


SYMBOL 


SENSOR 


Roll rate, c.g. 


'Peg 


Humphrey rate gyro. Model 01-0238 


Pitch rate, c.g. 


®cg 


Humphrey rate gyro. Model 01-0238 


Yaw rate, c.g. 


^cg 


Humphrey rate gyro. Model 01-0238 


Ro 1 1 rate, seat 


4’seat 


Minn.-Honeywel 1 rate gyro. Model 
JG 7005A4-9 


Pitch rate, seat 


®seat 


Minn.-Honeywel 1 rate gyro. Model 
JG 7005A4-9 


Yaw rate, seat 


'I'seat 


Mi nn.-Honeywel 1 rate gyro. Model 
JG 7005A4-9 


Pitch and roll rate, 
vertical gyro 


®VG 


Minn.-Honeywel 1 rate gyro. Model 
JG 7044A4 


Lateral acceleration 


ay 


Donner linear accelerometer. 
Model 4310 


Normal acceleration 


az 


Donner linear accelerometer. 
Mode 1 4310 


Rudder deflection 


Sr 


2 K Q, modified, spring loaded, 
selsyn potentiometer 


Left ai leron def 1 . 


Sa 


IRC 2 K Q dual, spring loaded, 
potentiometer 


Left spoiler defl. 


Ss 


Bournes 4 i nch 1 i near 
potentiometer 


StabI lator defl. 


Ss 


2 K Cl, modified, spring loaded, 
selsyn potentiometer 


Rudder pedal defl. 


Sped 


spring loaded, modified selsyn 
potent i ometer 



Lateral stick defl. 6gt}ck lat spring loaded, modified selsyn 

* potentiometer 



TABLE I 
(Conti nued) 



QUANTITY SYICOL 

Longitudinal stick 6 stick. 

Left and right 

throttle position , 6^ 

R 

Airspeed Vq 



■SENSOR 

long spring loaded, modified selsyn 
potentiometer 



spring loaded, modified selsyn 
potentiometers 



Altitude 



h 



W. 0. Leonard Model 50-2100 
W. 0. Leonard Wbdel 50-2000-9 
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TABLE II 

RECORDED FLIGHT DATA QUANTITIES, 
TELEMETERING CHANNELS UTILIZED 



QUANTITY 


CHANNELS UTILIZED 


Roll rate, eg 


3 


Pitch rate, eg 


3 


Yaw rate, eg 


3 


Roll rate, seat. 


3 


Pitch rate, seat. 


3 


Yaw rate, seat. 

Pitch rate, vertical gyro 
Roll rate, vertical gyro 


3 

1 

1 


Lateral acceleration 


2 


Normal acceleration 


2 


Rudder deflection 


2 


Left aileron deflection 


2 


Left spoiler deflection 


2 


Stabi lator deflection 
Rudder pedal deflection 
Lateral stick deflection 
Longitudinal stick deflection 
Left throttle 
Right throttle 
Airspeed 


2 

1 

1 

1 

1 

1 

1 



A Iti tude 
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TABLE III 

WIND CONDITIONS DURING CARRIER TEST APPROACHES 
(Headwind and Crosswind Components) 



RELATIVE RELATIVE 

APPROACH NUIVBER HEADWIND (KNOTS) CROSSWiND (KNOTS) 



1 


36 


0 


2 


36 


0 


3 


36 


0 


4 


36 


0 


5 


36 


0 


6 


36 


0 


7 


36 


0 


8 


19 


0 


9 


19 


0 


10 


19 


0 


II 


29 


0 


12 


26 


0 


13 


37 


0 


14 


31 


0 


15 


28 


0 


16 


26 


0 


17 


25 


0 


18 


25 


0 


19 


25 


0 


20 


23 


0 


21 


23 


0 


22 


23 


0 


23 


23 


0 


24 


26 


0 


25 


24 


7 1/2 right 


26 


28 


5 right 


27 


24 


5 r i ght 


28 


27 


5 right 


29 


35 


0 


30 


34 


6 right 


31 


34 


5 right 


32 


34 


9 r i ght 


33 


30 


5 r i ght 
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table IV 

RMS VALUES OF TURBULENCE QUANTITIES 



APPROACH 


“9 


9g 


Pg 


«Pg 




(doarees) 


(dea/sfic) 


(decrees) 


(dec/ sec) 


(dec/ sec I 


1 


2.21 


3.06 


1.94 


2.88 


3.39 


2 


2.21 


2.78 


1.69 


- 


- 


3 


2.39 


2.89 


1.47 


2.77 


4. 18 


4 


2.21 


3.41 


- 


2.89 


3.88 


5 


2.16 


2.72 


1.75 


3.19 


3.55 


6 


2. 16 


2.77 


1.66 


3.27 


3.75 


7 

8 
9 

10 
1 1 
12 
13 


2.16 


3.17 

DATA 


1.79 

UNUSABLE 


3.27 


3.73 


14 


1.64 


- 


1.86 


2.70 


3.76 


15 


- 


3.60 


2.03 


2.71 


3.66 


16 


1.66 


- 


- 


3.34 


- 


17 


1.63 


2.69 


- 


- 


- 


18 


1.53 


3.17 


1.97 


- 


3.11 


19 


1.43 


2.99 


1.97 


3.03 


3.66 


20 


- 


- 


1.56 


3.34 


4.25 


21 


- 


2.97 


- 


2.93 


3.34 


22 


2.32 


2.97 


1.44 


3.33 


4.24 


23 


- 


3.17 


- 


3.22 


3.66 


24 


2.31 


3.14 


1.61 


3.51 


3.89 


25 


2.41 


2.89 


1.57 


2.95 


3.48 


26 


2.09 


3.34 


- 


3.04 


3.78 


27 


- 


3.16 


1.77 


- 


3.21 


28 


- 


- 


- 


- 


- 


29 


2.51 


- 


1.67 


2.84 


- 


30 


2.23 


3.17 


1.98 


3.14 


3.76 


31 


2.09 


3.27 


1.91 


2.86 


3.89 


32 


2.43 


3.17 


2. 10 


2.99 


3.67 


33 


2.20 


3.27 


2. 10 


3.04 


3.83 







I 




1 









(0 






1 











r4 




C CuRVE.^ 



I Sanborn f^tctyuHn^ f^e^^mapaptft I I 















MO'^VnOHdOa SIOHINO 



3UFFAL0. NtW YORK 



pRirviTrn in «< s a 






rr 












■rtr- 




I Sanborn f^exx^uilnq f^e^^rttapapeft I i 





(■# 











I Sanborn ffuxyuiinq HunaptY^ I 






« 

^^/iicorc/ 






N/iVtON P^eL/MfNj9Ay. 



-4 





t 

t 



>4 



,k 

o 



v> 



vl 

I 



k 




/ 




r- 



/=VV5 /f re^^&L. /^A' /^LOdj 



I 



I 




I I Sanborn Rtctytdinq 








i 



\ r \ 

lO ^ 

•''I . 

V!T^!o 
y Nl 

o <> 












liiiiiii 










GRAPHIC CONTROLS CORPORATION 






I buffalo, mew yOPK printed in u s a | 





r ON 











•I .1 



I 



1 



' i ' 




\- ~ ¥- i i l it < . < t I < » t i t ♦ i I I I i I i ♦ t -»-»■ i » 




PRINTED IN U S A 



I I ^ PAC o r^qTn g c h a RJT S 




( 




l i i i i r I 1 I i 4 t ^ >-t I l-i-t — t -l - t t- i t i i t t- t « ^ »— 



4 - . 




4-- t i I I I >-1-1 -1 I « > t -i^ t I I I I I » t - l — i- i i r -i-lt i i‘i i 





I BUFFALO. NEW YGRK PRiNreoiNusA | 





PRInOo in H j" g I I 

rr- iiU-4t‘* ' : • : • 



CORDING CHAR TSl I 




♦-t -t i l i i i - i I I I t -4 I t t— t— l-t t ■>■ I -»-l t i-4-i-4- 



■irf 




JG CHARTS I GRAPHIC CONTROLS CORPORATION 



I BUFFALO. NLW YORK 




# 




m 



GRAPHIC CONTROLS CORPORATION 



I 




♦-•- I I i j_f j - i -i ♦ I < I -I i i - t > t i I t i t * ♦ t t i I I I jt I 





I a I§ RE COR pi MG CHA'RT^ | GRAPHIC CONTROLS CORPORATION 







I 




)RP0RATI0N I BUFFALO. NEW YORK 




:ha rts' 1 I GRAPHIC COMTROLS CORPORA! iQN 






NFW v.)Rk 



PRiNTtD IN U S A 




GRAPHIC CONTROLS CORPORATION 



I BUFFALO NtW V. f^K 



PRINTt O IN U S A 



I 





1- 1 I i i t i - i- i - i - i i i i i i i t t-i i t i i- t - t - t i -t-l t > t t t 



i 





i 



I 



i 



® RECORDING CHARJ’^ | GRAPHIC CONTROLS CORPORATION 




P».itsi rt O IN U j A 



DN I buffalo, NtW YORK 







i » * t t 1 4 t k I t I I I < - > i * r i H -l-- t -*-1 i r - r-i 



I 




4-^ - 4 > i 4»i ♦ 4-fr i t i I t i > ■ < »- 1 i t »-4 t t t i • • 4 >i.4 irrin 



dRDING CHARTS ! | GRAPHIC CONTROLS CORPORATION 






I i^UFf^ALO. NLW YORK 



PRINTED IN U 5 ^ 



I 









■^2ae=t-& 




TTir rr*rTTT 



t k i i > I "t t I i »--i i l- t - i i r« » i * i -i t t-l i i i i i - 4 



I 





4 i > -l i i - l t l t >- t - 4 i !-»■ I >■ i 



I l:^ ^*^t._ORDLNG CHARTS] | GRAPHIC CONTROLS CORPORATION 




>ORATION I BUFFALO. NEW YORK PWlNTt O IN U S A I 





I BUFFALO, NEW YBJK 



TION 



PRINTED 





I BUFFALO, NbW YORK 



^OLS CORPORATION 







I 



t 






4, 



(g RECO_ROING J^H^RTS] | GRAPHIC CONTROLS CORPOflAl ION 



I Ml^^^ALO, NtW YORK 








PRINTED IN USA | . | ♦ g R E C O H OiN G C H A R 




s] I GRAPHIC CONTROLS CORPORATION | BUFFALO. NLW YORK ph.nuoinusa | 




I 



1 



J 




till! il H — I I I - I- I I » -I t ♦ l-t t-t t I II I lut l-l-li jrt l i . 








LQHARTS] I graphic controls corporation 



I RUFFALU Ntw YORK 



PHtNTtO IN u S A 




♦ I I I t I I t >-|- t t I I > t I I I » t < * I t > * « ■« » « « I » 



I 

! 

I 



I 






( 




^R ECOR OIWG CHARTS! | GRAPHIC CONTROLS CORPORATION 



BUFf ALC 









l -i III- I t I I I t I -I li lt i t I li > I I » » 




I BUFFALO. NEW YORK P«.NttD.NusA | 





I 






e 



I 




^ PE COR^WG C HARTS | GRAPHIC CONTROLS CORPORATION | BUFFALO. NEW YORK 







[^RECOROiNG CHARTS] | GnAPHIC CONTPQlS COHr’HHA ' :GN 



PUFPAlO fVJbW YOPk 



I 



o 



i 





(^RECORDING CHARTS] | GRAPHIC CONTROLS CORPORA,' 






t I I ! ■ > l- t I l - t I I I I t I I I t I I I t 



5QRAH0N I BUFFALO. NEW YORK 




M t I I t i •%-% t I I I I i t -t*-! I I » < < ♦ "< * i - t " » t- f i t-t 




iRPORATION I BUFFALO. NEW YORK 






d 




• ■* • 






I ^IECOROLNCj CHART'S] | GRAPHIC CONTROLS CORPORATION | BUFFALO. NEW > 





♦t I I I I I t I t -i I I t- t - l -l -l I • I I • I I -I I l l-t i til l . 









I 




■ 





IgR ECORDI NG ^ARTSl | GRAPHIC CONTROLS CORPORATION 



I RUFf-ALO. NIW< 



:I44_’4^ 



^ I I I I t t i I I I I I » l-»- • l - t t I t I t H I t--4 f ^ 




ITEO IN U S 



^ RECORDING CHART S] | GRAPHiCj 




iH - I I I I i I ^ t I r i I I t I »-•-! > -i i H > 




J 





PHIC CONTROLS CORPORATION | BUFFALO. NEW YORK printcdinusa | 




5 1-51 



^ RECdRDIN G CHARTSI | GRAPHIC CONTROLS CORPORATION 

/¥ 




-♦ I lil t i- i-i - f-i I t I I- I I I I- 1 I I I 1 I f i -f -n I I I t- t i % 






I 



1 



[^RECORDtNG CHARTSl | GRAPHIC CONTROLS CORPORATION 



I buffalo. Nt W YORK 




f-f -I -I- 1- 1 -f- t- i-i i-t-i-i I t I I • t t I i i i I I i- t - t i I i -i t L 




I 








NTROlS CORPORAT-UtN 



-ato! 




-T— n 

lati- 


rr 

L^.r 




ii, . I 




1 '-1 ■ 





4 1 >- t - t t t -e i i i -f 8 t t »-i ■ t - t-t ♦-» 



I BUFFALO. NEW YORK 




4 -i ri -i I I i I I t -t < I t -l- i t t < t I i i-4-i t » i t < l i- L 



i RECORDIIMG CH ARf^ | GRAPHIC CONTROLS CORPORATION 




t t t t t ri i" i i i i i 4 . 4n |_ 



BIrTSI I GRAPHIQ^Or^OLS CORPORATION 

Tiafc <xfT£^ 





C- I {^RECORDING CHA'RtS | u\PH< iONTROLS mRPi HATlON 




I 



^ CORPORATION 



I BUFFALO. NEW YORK 



PRJNTtO IN U S A 




1 




RECO RDiNG CHARTS] | GRAPHIC CONTROLS CORPORATION 







1 

: ., / 


1 1 
'1 '1 

1^' 


1 

r ■-'_ 






■1 


V 

- [ 


1 ' 
.1 


If 






i 



I -I I > » I I--I I I I -t I l- t-l I I I -A t- t 1 -1 I -I — t-i - 1— I - t — ♦— l - l - I 




[g R E C O R q » IM G ‘ C R T S ] | GRAPHIC CONTROLS 



CORPORATION 



I RUFFAlO NtW YlJRk 





CHARTSI I GRAPHIC CONTROLS CORPORATION 



I BUFFALO. NtW YORK 



PRINTC O IN U S A 




( 



^ I ^ ^ 65N54 I (gREC6R DIMG CI T arT^ | GRAPHIC CONTROI T, ^ 





I 





m 

I 




USA 




♦- I - t I I I I t I t I < I t I — f- t -f- t - i t I i t t I t i- t I t I -f-i - i- f 



I BUFFALO. NLW rOBK 

1 





! 

J 







i 









IgRECORmNG charts! I L.HAPHJC CUNlHOlb 




I * I I ■» I »-i I t t I I i i t I t I i -t- • I- - I <1 l i i » * t"-« » t 




1 




I BUFFALO. NEW YORK printed in u s a | 









I 

I 



H RE^CQRDING ch arts] I graphic controls corporation I buffalo, new YORK 




J 



D. NEW YORK 



USA 



t> .> > 



I 



gRiea 



I 





^ o s fK I I (§ RE CORDING CHARTS^ | GRAPHil. ^;0^gTR; 




Yl iMK 



^ R E C O R DiRJ G t H AJ^ T 6 1 ( GRAPHIC CONTROLS - IRPORA ^ )N 




I hljfau 




Nr v\ 



PWINTI O tN 







' 1 


i ■ : 

, . I 













» I 





i I t i . | t % I I t I I I f t I t t I l‘t m I i _ t 



R E Cd R Dj N G C M A R | GRAPHIC CONTROLS LUHPuRATlON | BlJf-FAnJ NLW YORK 




m 



« 




(§ R E 1 0 R a I W G C HARTS] | GMAPHI . CONTROLS 



PORAI iON 






<3urr» 








1 



i 

I 





— f 




I 




BUFFALO. NFW YORK 



PRINTED IN U S A 



I BUFFALO. NEW YORK panted in us a | 




I 




I I BUFFALO NEW YORK pR.NTtoiNusA , 




^?7/f^/sZ Sof'S/^ 



/at t/^r^cus" /‘f^^roscheS' , 




/hese /e>ut ^^radch-es" 
z--^/ fb )C-u>/rKJ^, 



'i9to )C-£o/r?c/^ 
ZWa-^s/ 
a/v 7 ^ 

IS ^f?o/s) 



£f ’/'0 
^/oyy^ 

Kt)o/x J^ 




I 



thes037 

The acquisition, reductin, and analysis 




3 2768 001 97477 7 
DUDLEY KNOX LIBRARY 



